Temporal variability of the telluric sodium layer 
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The temporal variability of the telluric sodium layer is investigated by analyzing 28 nights of data obtained 
with the Colorado State University LIDAR experiment. The mean height power spectrum of the sodium layer 
was found to be well fit by a power law over the observed range of frequencies, 10 /iHz to 4 mHz. The best 
fitting power law was found be be W^u 01 , with a = -1.79 ±0.02 and (3 = 1.12 ±0.40. Applications to wavefront 
sensing require knowledge of the behavior of the sodium layer at kHz frequencies. Direct measurements at 
these frequencies do not exist. Extrapolation from low-frequency behavior to high frequencies suggests that this 
variability may be a significant source of error for laser-guide-star adaptive optics on large-aperture telescopes. 
© 2008 Optical Society of America 
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1. Introduction 

The advances made in adaptive optics (AO) over the last 
decade have enabled many ground-based optical astro- 
nomical projects that were never before possible. 1 Most 
of these studies have been performed with natural guide 
star AO systems (NGSAO), which require a bright star 
near to the desired target. While there are many exam- 
ples of strikingly successful studies done with NGSAO, 
the fraction of the sky accessible to them is rather small. 
In order to access the whole sky, one requires the abil- 
ity to sample the atmospheric distortion anywhere, not 
just toward the direction of a bright star. The telluric 
sodium layer is a layer of sodium atoms that is typi- 
cally present at an altitude of ~ 90 km over an extent of 
~ 15 km. By illuminating the layer with a laser tuned to 
the wavelength of the Di transition (589 nm), one may 
excite emission from these sodium atoms. The resulting 
bright spot acts as an artificial guide star, and enables 
the operation of laser guide star AO (LGSAO) systems. 

The development of a LGSAO system is of particu- 
lar importance to the Thirty-Meter telescope (TMT) . In 
order to take full advantage of the potential spatial res- 
olution of the large aperture, it is imperative for many 
of the TMT's instruments to be fed by the proposed 
LGSAO system, the Narrow Field InfraRcd Adaptive 
Optic System (NFIRAOS). The development of NFI- 
RAOS presents many technical challenges. As discussed 
in O'Sullivan et al. (2000), 2 one of these challenges — 
the refocusing of the telescope on the moving laser guide 



star spot — is due to the variability of the sodium layer. 
Though the physical mechanism is not entirely under- 
stood, atmospheric gravity wave propagating upward in 
the mesopause region is thought to be the most likely 
explanation. A typical Shack-Hartmant wavefront sen- 
sor will determine the centroid of a illuminated spot on 
a CCD, and will therefore respond to the instantaneous 
mean height of the sodium layer. A change in the mean 
height, Aft, of the sodium layer produces a RMS wave- 
front error, averaged over the aperture, 
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where D is the diameter of the telescope aperture, and ft 
is the instantaneous density- weighted mean height of the 
sodium layer. 3 Note that a is proportional to the square 
of the diameter of the telescope. Thus, for the TMT, we 
expect the effect to be an order of magnitude larger than 
for the largest optical telescopes presently in operation. 
For the TMT, a = 8 nm for every meter change of the 
mean height. 

While an AO system can use laser and natural guide 
stars to correct for focus error caused by the atmo- 
sphere and variations in the height of the sodium layer, 
it can not do so instantaneously. Therefore, the correc- 
tion applied by the AO system will be imperfect, and 
a residual wavefront error (RWFE) will remain. While 
the mean height of the sodium layer is known to be vari- 
able on the order of several kilometers over timescales of 
minutes, LGSAO systems sample the atmosphere many 
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Fig. 1. shows the time series of the centroid height for 
December 24 th . Channel #1 is shown as a solid line, 
while channel #2 is shown as a dotted line. 
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Fig. 2. shows the mean height power spectrum and 
power-law fits for December 24 th , with a indicating the 
slope, and (3 indicating the offset. 



times per second. Therefore, it is necessary to know how 
the mean height of the sodium layer will vary on millisec- 
ond timcscales in order to estimate the RWFE. In order 
to assess the variability of the sodium layer on short time 
scales we examined data taken by the Colorado State 
University (CSU) Light Detection and Ranging (lidar) 
experiment over 89 nights in 2003. We determined the 
typical RWFE expected by multiplying the mean height 
power spectrum of the sodium layer by a model rejection 
transfer function determined specifically for the TMT. 

2. Technique 

A. Reduction and extraction 

The CSU lidar samples the sodium layer by projecting 
laser pulses upwards, cyclically tuned to three prese- 
lected frequencies within the Doppler-broadened band- 
width of the sodium D2 transition at 589 nm. The inten- 
sity of the return signal in three channels (one for each 
frequency), binned as a function of distance from the 
telescope can be processed for measuring atmospheric 
temperatures and horizontal wind between 80 and 105 
km. 4 The experiment employs two telescopes that can 
be independently pointed, and can therefore measure 
the Na-layer in two directions simultaneously, typically 
zenith, due north at an altitude of 60°, or due east at 
an altitude of 60°. In order to increase the signal-to- 
background ratio, the signal is integrated for a duration 
of 120 seconds. Due to readout time, measurements are 
typically spaced on 125 second intervals, resulting in a 
Nyquist frequency of 4 mHz. 

The signal-to-background ratio was calculated in the 
following manner: The atmosphere was divided into 
three regimes. The signal regime extended from 75 km 
to 105 km. Although the sodium layer only extends over 
approximately half this altitude range at any one time, 
the large range was necessary to ensure that all the sig- 
nal was contained in the range regardless of the mean 
height. The two background regimes extended from 60 



to 72 km and 108 to 120 km respectively. By integrating 
over the aforementioned regimes a total signal, S to tah 
and total background, B to tah were calculated. The re- 
quirements for acceptable data were S total 

> 10 4 pho- 
tons and B to tai < 10 3 photons. 

The raw files were processed to yield time, mean height 
(corrected for zenith angle of the laser, and the altitude 
of the CSU lidar experiment), and a flag indicating if the 
signal-to-background ratio of an individual measurement 
was acceptable. The time series for December 24 th is 
shown in figure^ Channel #1 and channel #2 arc shown 
as solid and dotted lines respectively. 

B. Transformation and Fitting 

We computed the mean height power spectrum for each 
night individually. Over the observed range of frequen- 
cies, from /iHz to mHz, the transformed data were well 
fit by a power law. Many of the nights were split into 
several runs of acceptable data quality separated by sec- 
tions of poor-quality data. The extracted power spec- 
trum is a convolution of the Fourier transforms of both 
the data and the window function. A window function 
with a broad Fourier transform will yield a power-law fit 
with an artificially shallow slope, which will in turn pre- 
dict an inflated value for the RWFE at high frequencies. 
In order to derive the correct slope of the power-law fit 
we used only the longest continuous series of measure- 
ments from a given night. The power spectrum for the 
night of December 24 th is shown in figure [21 Channel #1 
and channel #2 are shown as solid and dotted lines re- 
spectively. The slope of the power-law fit for the various 
channels are a and the offsets are /3. 

C. Determination of residual wavefront error 

Because the distance to the sodium layer is not know 
a priori, natural guide stars are required to keep the 
telescope focused on science targets. In order for a suffi- 
cient number of guide stars to be present to make a large 
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Fig. 3. shows the cumulative distribution for the cal- 
culated RWFE for both channels of the 28 nights with 
more than 7 hours of data with acceptable data quality. 

fraction of the sky available, the NGSAO system must 
be designed to work with relatively faint stars. The low 
fluxes from these stars necessitate a slow read out rate. 
Therefore, NGSAO can only provide low-frequency infor- 
mation regarding focus changes. A laser guide star can 
be designed to emit relatively high flux, and can there- 
fore be read out quickly, providing high-frequency focus 
information. Herriot et al. 3 combine the high-frequency 
focus information from the laser guide star and the low- 
frequency focus information from the natural guide star 
to obtain focus information over a wide range of fre- 
quencies. Due to limited bandwidth, the aforementioned 
correction is not perfect, but may be characterized by 
the rejection transfer function. The residual wavefront 
variance (RWFE) 2 is proportional to the integral over 
the appropriate range of frequencies (10 -4 to 10 4 Hz) of 
the product of the mean height power spectrum and the 
square of the AO system rejection transfer function, 
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The residual wavefront error was calculated for each 
night with a rejection transfer function designed for the 
TMT. The cumulative distribution for all nights with 
more than 7 hours of continuous acceptable data is 
shown in figure |3J We found that nights with less than 7 
hours of continuous data had a shallower slope than the 
datasets with longer continuous streams. This is due to 
the aforementioned effect of the convolution of a narrow 
window function with the true power spectrum. 

3. Results and Conclusions 

The data for each night with at least one hour of accept- 
able quality (76 of 89) were analyzed. It was found that 
nights with sections of usable data that were less than 
7 hours (48 nights) had a shallower slope with respect 
to those nights with longer sections of usable data (28 
nights). Using just the longer nights we find the best- 



fit power law has the parameters a = —1.79 ± 0.02 and 
P = 1.12 ± 0.40. We found that RWFE for those nights 
has a median value of 20 nm and a 90 th percentile value 
of 34 nm. This effect is an important contribution to the 
overall error budget of NFIRAOS, and more generally, 
will likely be for most, if not all, of the LGSAO systems 
designed for the next generation of large ground-based 
telescopes. 

This analysis was limited by the relatively long in- 
tegration times of the CSU lidar experiment, which in 
turn set a low Nyquist frequency. The variability of the 
sodium layer height in the frequency range observed is 
mainly due to the upward passing of atmospheric waves. 
Atmospheric turbulence at the sodium layer height can 
perturb the mean height of the sodium layer on time 
scales of seconds or shorter. However, changes caused 
by atmospheric turbulence are likely to be small (cm or 
less), unlike the kilometer-scale changes reported here. 
The assumption that the mean height power spectrum 
can be fit by a single unbroken power law from /iHz to 
kHz is perhaps naive, but reasonable. The actual high- 
frequency behavior of the mean sodium height power 
spectrum, of course, can only be determined by further 
experimentation. The temporal resolution could be im- 
proved by designing a lidar experiment with either a 
larger collecting telescope, or by sacrificing spatial reso- 
lution. 
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